Azimuthal correlations of electrons from heavy-flavor decay with hadrons in p-\-p and 

Au+Au collisions at y / s/v r /v =200 GeV 
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Measurements of electrons from the decay of open-heavy-flavor mesons have shown that the yields 
are suppressed in Au+Au collisions compared to expectations from binary-scaled p+p collisions. 
These measurements indicate that charm and bottom quarks interact with the hot-dense matter 
produced in heavy-ion collisions much more than expected. Here we extend these studies to two- 
particle correlations where one particle is an electron from the decay of a heavy-flavor meson and 
the other is a charged hadron from either the decay of the heavy meson or from jet fragmentation. 
These measurements provide more detailed information about the interactions between heavy quarks 
and the matter, such as whether the modification of the away-side-jet shape seen in hadron-hadron 
correlations is present when the trigger particle is from heavy-meson decay and whether the overall 
level of away-side-jet suppression is consistent. We statistically subtract correlations of electrons 
arising from background sources from the inclusive electron-hadron correlations and obtain two- 
particle azimuthal correlations at ^snn =200 GeV between electrons from heavy-flavor decay with 
charged hadrons in p+p and also first results in Au+Au collisions. We find the away-side-jet shape 
and yield to be modified in Au+Au collisions compared to p+p collisions. 

PACS numbers: 25.75.Bh, 25.75.Gz 



I. INTRODUCTION 

Experiments at the Relativistic Heavy Ion Collider at 
Brookhaven National Laboratory have produced a hot 
dense partonic matter [3, Q]. Results from high-pr tt° 
production indicate that fast partons moving through 
the matter lose a substantial amount of energy through 
interactions 0, H|. This energy loss was expected to 
be reduced for heavy charm and bottom quarks due 
to the dead cone effect which suppresses gluon radia- 
tion However, electrons from the semileptonic de- 
cay of D and B mesons are seen to be suppressed at 
nearly the same level as 7r°s out to the highest mea- 
sured pr, ~10 GeV/c @. This challenges the picture of 
gluon radiation as the dominant means of parton energy 
loss. Various alternative scenarios including collisional 
energy loss Q, in-medium formation and dissociation of 
the heavy meson or an increase in the fraction of 
heavy quarks carried by baryons [tl [Io| have been pro- 
posed to account for the large suppression. 

Single particle yield measurements provide information 
on the overall deviation of particle production from p+p 
expectations. However, the observed high-p^ spectra are 
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thought to be dominated by particles that have lost less 
than the average amount of energy, either due to a short 
path length through the matter or by a fluctuation. In 
order to get more detailed information about the interac- 
tions between the particles and the matter two-particle 
azimuthal correlations have been extensively used. In 
p+p collisions these correlations are characterized by two 
back-to-back jet peaks At small azimuthal angular 
difference, A<f>, particles are from the fragmentation of 
the same jet; at A(f> ss tt particles are from the fragmen- 
tation of partons in the opposing jet. 

In heavy-ion collisions, these correlations can provide 
information about the pattern of energy loss for the back- 
to-back dijet system as well as other interactions be- 
tween the fast partons and the medium. Measurements 
of hadrons associated with a high-p^ hadron have shown 
the away-side correlations from back to back dijets to be 
significantly suppressed [l2l - fl5l |. Measurements of the 
correlations of electrons from heavy-flavor decay with 
other hadrons in the event can also provide insight into 
heavy-flavor energy loss and how this compares to tt° 
and direct photon triggered correlations where the mod- 
ifications could be different due to the different partons 
probing the matter. This is crucial for building a quan- 
titative understanding of the nature of the interactions 
between hard partons and the produced hot matter. 

In addition, a strong broadening and double peak 
(shoulder) structure of away-side correlations at moder- 
ate pt has been observed [16l |, Many theoretical ideas 
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have been proposed to explain this modification, includ- 
ing Cerenkov gluon radiation [I?], EH, large angle gluon 
radiation [l^, [20I and Mach shock- waves [2l|. Measure- 
ments of the shoulder structure with particles from the 
fragmentation of heavy quarks, especially bottom, are 
interesting because at moderate momenta the quark ve- 
locity will be much smaller than the speed of light, in 
contrast to light quarks where v «c at all jet momenta. 
In a Mach shock-wave scenario the cone angle of the dou- 
ble peak structure away from n, {8m), is related to the 
speed of the parton by: cos 9 M — ^ where eg is the speed 
of sound in the matter and v is the speed of the parton as 
it propagates through the matter. It has been proposed 
that double-peaked correlations that do not obey Mach's 
Law could favor strongly-coupled AdS/CFT string drag 
scenarios [22J , or transverse flow [23] . An alternative ex- 
planation based on geometrical fluctuations in the initial 
state leading to triangular flow has also recently been 
proposed [24| - |26j . 

Measurements in p+p collisions are a necessary base- 
line to heavy ion measurements, particularly for heavy- 
flavor triggered correlations. At leading order, several 
sub-processes contribute to charm production leading to 
a midrapidity D meson. For px <10 GeV/c, one lead- 
ing order calculation shows that ~20% of the time the 
charm quark leading to the D, which decays semileptoni- 
cally into the trigger electron, is balanced by an opposing 
(A(j) a; 7r) c quark |27| . The rest of the contribution is 
from processes such as eg — ¥ eg or cq(q) — > cq(q) where 
the c is not balanced by a midrapidity, high-p^ c. Next 
to leading order effects are known to be large in heavy 
quark production. The powheg Monte Carlo calcual- 
tion |28j |. which includes 2— >3 processes, also shows sub- 
stantial contributions to the away-side correlations from 
gluons. Thus, in order to measure cc or bb correlations 
one should identify the heavy quark in both the trigger 
and away jets. For the present purposes, this means it is 
not possible to identify the jet opposing the electron from 
heavy-flavor decay unambiguously as also from heavy- 
quark fragmentation. This complicates the interpreta- 
tion of the present measurements. However, the com- 
parison of heavy- and light-flavor triggered correlations 
can still provide a crucial step toward understanding fast 
parton propagation through the matter. 

In the present work heavy-flavor electrons are those 
from the decay of both D and B mesons. The relative 
contribution of electrons from bottom to the total heavy- 
flavor electron yield changes with the pr of the electron 
and has been measured in p+p collisions to be from wlO— 
50% for 1.0< p T <6.0 GeV/c (HHj. Fixed-order plus 
next-to-leading-log (FONLL) calculations from Ref. [3l[ 
agree well with the measured bottom contribution. 

We present first results of the azimuthal correlations 
between electrons from the decay of heavy-flavor mesons 
with charged hadrons in Au+Au and p+p collisions. 
We statistically subtract correlations from electrons due 
to background electron sources (Dalitz decays, photon 
conversions and quarkonia) from the measured inclusive 



electron-hadron correlations. 

The paper is organized as follows. In Section|TT]we out- 
line the analysis procedure used; in Section [TTT] we show 
the results in p+p and Au+Au collisions; and in Sec- 
tion |IV] we conclude and discuss the prospects for future 
measurements. 



II. ANALYSIS METHOD 
A. Experimental Setup 

These results are based on 1.1 billion level- 1 trig- 
gered p+p events sampling 8.0 pb~ x taken during the 
2006 RHIC running period and 2.6 billion minimum bias 
Au+Au events, corresponding to 0.41 nb~ 1 taken dur- 
ing the 2007 RHIC running period. The events were 
triggered by a hit in each of two Beam-Beam Counters 
(BBC) at 3.1< \i]\ <3.9 and the interaction is required 
to be within 25 cm of the center of the interaction re- 
gion. The p+p level- 1 triggered sample also required an 
energy deposit of approximately 1.4 GeV in an overlap- 
ping tile of 4x4 EMCal towers in coincidence with the 
BBC trigger. EMCal towers are A(j) x Arj w 0.01 x 0.01. 
In Au+Au collisions the event centrality is measure d by 
the charge seen in the Beam-Beam Counters (BBC) [32j |. 

The charged particle tracks and photons used in this 
analysis are measured in the PHENIX central arm spec- 
trometers. Electrons are measured between 1.5 and 
4.5GeV/c and charged hadrons are measured between 
0.5 and 4.5GeV/c. PHENIX has two such spectrome- 
ters, East and West Arms, each covering tt/2 rad in az- 
imuth and <0.35. This analysis uses in each arm a 
drift chamber (DC), two layers of pad chambers (PCI 
and PC3), a Ring Imaging Cerenkov Detector (RICH), 
and an electromagnetic calorimeter (EMCal). Charged 
particles (both electrons and hadrons) are reconstructed 
in the DC and PCI. Electron identification is done by 
requiring two (three) associated hits in the RICH for p+p 
(Au+Au), a shower shape cut in the EMCal and an E/p 
cut, where E is the energy of the cluster in the EMCal 
and p is the track momentum determined by the DC. 
Electron candidates are required to have a matching hit 
in the EMCal within 3er (2a) in p+p (Au+Au). Cuts on 
the RICH ring center and shape are also included. The 
hadron contamination remaining is less than 1% in p+p 
collisions and less than 3% in Au+Au collisions. Hadrons 
are identified by a RICH veto and a confirming hit in the 
PC3. Photons are identified by a shower shape cut in the 
EMCal and a veto in the PC3 to reject charged tracks. 
Hadron contamination in the photon sample is less than 
4%. Cuts on electron-hadron and photon- hadron pairs 
are also used to equalize the pair acceptance between 
real and mixed pairs and remove pairs that share hits in 
the various detector subsystems. 

In the Au+Au running period the Hadron Blind De- 
tector (HBD) [33d was installed for a commissioning run 
between the beam collision vertex position and the cen- 
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tral arms. Photon conversions in the detector material 
were an additional source of background electrons in the 
inclusive electron sample. The HBD in front of the West 
Arm was absent for a substantial portion of the running 
period. We make the additional requirement that elec- 
trons from the 2007 running period are reconstructed in 
the West Arm and select events only from the running 
period where the HBD in front of the West Arm was 
removed in order to reduce the number of photon con- 
versions. 



B. Background Subtraction 

In two particle correlations there is a combinatorial 
background due to pairs where the particles are uncorre- 
cted except by event- wise correlations, such as central- 
ity and the reaction plane in Au+Au or the underlying 
event in p+p collisions. This background is very large in 
central Au+Au collisions and much smaller in p+p col- 
lisions. In Au+Au collisions the background is removed 
by the Absolute Background Subtraction technique 34 1 
while in p+p collisions by the Zero Yield at Minimum 
(ZYAM) method [35| (with the uncertainties determined 
as in Ref. [34]) is used to subtract the A0 independent 
underlying event. In the Absolute Background Subtrac- 
tion method the combinatorial background yield is de- 
termined from the centrality dependence of the single 
particle yields. The principal advantages of this method 
for this analysis are that the background uncertainty is 
not subject to statistical fluctuations caused by the small 
number of electron-hadron pairs and no assumption is 
made about the shape of the pair Acf> distribution. 

In p+p collisions the underlying event independent of 
Acf>. However, in Au+Au collisions there are additional 
correlations due to elliptic flow, V2- These correlations 
do not affect the magnitude of the combinatorial back- 
ground, but do affect the azimuthal correlation of par- 
ticles. The t>2 values used in this analysis are measured 
using the reaction plane as determined from the PHENIX 
Reaction Plane Detector and the same particle cuts used 
in the correlation analysis. 

Results are reported as conditional yields of hadrons 
associated with trigger electrons after combinatorial 
background subtraction, Y e _^(A</)): 



y e _ h (A</») = 



1 dN„ 



N e e h dA(j) 



<-fl(l + 2wfu£cos(2A0)) (1) 



where N e is the total number of observed trigger parti- 
cles and €h is the reconstruction efficiency for the associ- 
ated hadrons as determined by a GEANT based Monte 
Carlo simulation and embedding single particles into 
real events. rf ^g^ as is the measured trigger-associated 
particle A(j) distribution, which has been corrected for 
nonuniform two particle A(f> acceptance by using mixed 
events 34] . B is determined by the background subtrac- 
tion methods described above. 



Removal of Non-Open Heavy-Flavor 
Electron-Hadron Correlations 



Studies of electrons from open heavy-flavor decay are 
complicated by the background of electrons from light 
meson decay, photon conversions and, at higher px, 
quarkonia and Drell-Yan. In this analysis we statisti- 
cally subtract the correlations from these background 
sources using a method similar to that used to mea- 
sure direct photon-hadron correlations [3t| . The yield 
of inclusive electron-hadron pairs per electron trigger, 
the conditional yield Y ein< ._h{PT,e,PT,h, A<j)), is measured. 
This is a weighted average of the conditional yield of 
hadrons associated with electrons from heavy-flavor de- 
cay and the conditional yield of hadrons associated with 
electrons from the background sources (the dependencies 
on pt and A(j) are omitted hereafter for simplicity) : 



6 bkg bkg 



N K 



(2) 



where N eilF {N ebh ) is the number of electrons from 
heavy- flavor decay (background sources). Y enF —h can 
then be written as: 



(fl HF + l)Y eznc - h - Y ebkg 

i?HF 



(3) 



where -Rhf 



N, 



The two quantities to be deter- 
i?HF can be determined 



mined are -Rhf and Y Bbk _ 
by comparing the measured electron yields to the known 
sources of background electrons. The Y &bkg -h determi- 
nation is based on measured and simulated azimuthal 
correlations of the sources of background electrons and 
is described in detail below. 

The i?HF value giving the composition of the elec- 
tron triggers into heavy-flavor and background sources is 
taken from the published PHENIX measurements 0, H3] 
for the electron px bins used in this analysis. Based on 
simulations the -Rhf has been decreased to account for 
extra air conversions due to the removal, in both the p+p 
and Au+Au data samples, of the helium bag, which was 
installed during the data taking periods of Refs. d, H3|- 
The removal of the He bag added 0.65% of a radiation 
length to the material in front of the tracking system. 
However, the reconstruction efficiency of the electrons 
from these air conversions decreases with the distance 
from the interaction point. This was simulated using a 
GEANT-based description of the PHENIX detector and 
the electron-identification cuts described above. Addi- 
tionally, electrons due to quarkonia decay are included 
in the background sample |38|, further reducing i?HF- 
The reduction due to quarkonia is determined by J/^t 
measurements and simulations of the decay [38l - l4fj ] . The 
change to the background yield due to quarkonia is 3% at 
1.5< p T <2.0 GeV/c and 38% at 4.0< p T <4.5 GeV/c in 
p+p collisions (at moderate and high pt the heavy-flavor 
signal is larger than the background, so the change to 
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TABLE I: Ratio of heavy-flavor electrons to background elec- 
trons in p+p along with the systematic uncertainty. 



ship can be written as: 



p T (GeV/c) 



Re 



1.5-2.0 
2.0-2.5 
2.5-3.0 
3.0-3.5 
3.5-4.0 
4.0-4.5 



0.66 ± 0.13 
0.86 ± 0.14 
1.09 ± 0.14 
1.31 ± 0.16 
1.49 ± 0.17 
1.60 ± 0.17 



TABLE II: Ratio of heavy-flavor electrons to background elec- 
trons in Au+Au along with the systematic uncertainty. 



p T (GeV/c) 



Rhf 



1.5-2.0 
2.0-2.5 
2.5-3.0 
3.0-3.5 
3.5-4.0 
4.0-4.5 



0.94 ± 0.21 
1.14 ± 0.25 
1.29 ± 0.30 
1.38 ± 0.34 
1.43 ± 0.36 
1.38 ± 0.36 



the heavy- flavor electron spectra is much smaller 1381]). 
At all but the highest px used here the change is only 
due to electrons from J /iff. Electrons from T become 
important at pt >4.0 GeV/c and electrons from Drell- 
Yan are negligible at all transverse momenta used in this 
analysis. The -Rhf values used in this analysis are shown 
in Table U (p+p) and Table [ED (Au+Au). 

The remaining unknown in Eq. [3] is the azimuthal 
correlations of the background electrons with hadrons, 
Y ebk h- These pairs can be divided into two classes: 
those from photonic sources, electrons from the decay 
of light mesons and photon conversions in the detector 
material, and those from quarkonia (electrons from light 
vector meson decay are a small contribution and are ne- 
glected). These correlations are determined from inclu- 
sive photon-hadron correlations. Inclusive photons, like 
photonic electrons, are largely from 7r° and r\ decay at 
Pt <5 GeV/c, and pythia (version 6.421) [41j simula- 
tions of the correlations between electrons from J de- 
cay with hadrons. The fraction of hadrons misidentified 
as electrons is small, as discussed above, and those corre- 
lations are not separately subtracted. At pr, e <2 GeV/c 
there is a small (~3%) contribution from electrons from 
the semileptonic decay of kaons whose correlations are 
also neglected. 

a. Photonic Electron Correlations Photonic elec- 
tron sources include Dalitz decays and photon conver- 
sions where the photons are from light meson decay. To 
determine these correlations, Ye. ot -h, we measure inclu- 
sive photon-hadron correlations. Inclusive photons are 
dominantly also from light meson decay. However, the 
parent meson pt distributions need not be the same for 
photonic electrons and inclusive photons. The relation- 



Y,. 



-h(PT,i) 



E 
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Wi(PT,j) Y-n„ a -h(pT,j) (4) 



where each i (j) represents a 0.5 GeV/c bin in electron 
(photon) pt- The weight coefficients Wi(pT,j) are deter- 
mined via simulation and are used to transform the in- 
clusive photon-hadron correlations into expectations for 
photonic electron-hadron correlations. 

Two methods are used to determine the Wi (ptj ) ■ The 
first method treats the electrons as coming from pho- 
ton conversions in the detector material and the second 
treats the electrons as coming from Dalitz decays. True 
photonic electrons come from both sources, however both 
methods give very similar uii(pTj) values. 

In the first method the measured single inclusive pho- 
ton spectrum is input into a GEANT based simulation 
of the PHENIX detector. The same electron identifica- 
tion cuts as in the real data analysis are then applied 
to reconstructed conversion electrons and the relation- 
ship between the input photon pt and the reconstructed 
conversion electron pt determines w. 

In the second method, the n° spectrum from Ref. [42j 
(for p+p collisions) or Ref. [|| (for Au+Au collisions) is 
taken as input to a Monte Carlo simulation that decays 
the 7r°s via Dalitz decay. The relationship between the 
intermediate low mass virtual photon and the resulting 
decay electron are used in the same manner as in the first 
method to determine w. Since the mass of the virtual 
photon is small the difference in the pt distribution be- 
tween real photons and the virtual photons in the Dalitz 
decay is negligible. 
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FIG. 1: (color online) Weights (w) for electrons with 1.5< 
pr <2.0 GeV/c from the method using Dalitz decays (his- 
togram) and photon conversions (solid points) as a function 
of the photon pt- 

Figure [T] compares the w from the two methods for 
a single electron pt selection and shows that the differ- 
ences between the two methods are small. The maximum 
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FIG. 2: (color online) e^^g — h (top panel) and 6hf — ft 
(bottom panel) conditional yields for p+p collisions with 
1.5< p T ,e <2.0 GeV/c and 1.5< p T ,h <2.0 GeV/c for the 
two methods of constructing e^ a — h conditional yields: the 
conversion method (solid circles) and the Dalitz decay method 
(black squares). Points have been offset slightly for clarity. 



deviation in the resulting enF-h conditional yields is at 
small A<f>, where the difference between the two meth- 
ods is 0.008 (0.006) in p+p (Au+Au). Both inclusive 
photons and electrons from Dalitz decays are largely the 
result of 7T° decay. The ir° spectrum falls steeply with px 
and thus the measured photonic electrons at a given p T 
are dominated by those carrying a large fraction of the 
7T° pt regardless of whether the electron comes from a 
conversion or a Dalitz decay. This argument holds for all 
heavier mesons that contribute to the photonic electron 
sample, which explains the small difference between the 
conversion method (including all mesons that decay into 
photons) and the Dalitz decay method (including only 
7T° decay). 

Figure [2] shows the difference in the ebkg — h and enp — h 
correlations for the two methods. The difference is small 
compared to the statistical uncertainty and is included 
in the systematic uncertainty. An additional small sys- 
tematic uncertainty is included from the statistical un- 
certainty of the simulations. 

b. Correlations of Electrons from Quarkonia Decay 
The azimuthal correlations between J/'Fs and hadrons 
have not yet been measured at these momenta, so 
PYTHIA [4l| is used to simulate the correlations between 
the electrons from J /iff decay with charged hadrons. For 
both the Au+Au and p+p measurements the default J/iff 
production within pythia is used. For p+p, the sys- 
tematic uncertainty is taken as the maximal deviation 
from the default production when varying the production 
mechanism between color singlet (pythia ISUB=421) 
and color octet (pythia ISUB=422) states. Figure [3] 
shows the correlations of electrons from J /iff decay and 
hadrons after ZYAM background subtraction for an ex- 
ample pt selection. For Au+Au, the situation is more 
uncertain as a substantial fraction of the J/'Fs could be 
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FIG. 3: (color online) pythia ej/y — h correlations after 
ZYAM subtraction for electrons with 3.0< pr <3.5 GeV/c 
and hadrons with 0.5< pr <1.0 GeV/c. The central values 
are from the default pythia J /iff production setting and the 
lines show the systematic uncertainty set by the magnitude 
of the maximal deviation between the default setting at the 
color singlet and color octet production settings. 



coming from rccombining c and c quarks [43l |44| . In 
this case, the azimuthal correlations could potentially be 
strongly reduced; the systematic uncertainty is taken to 
extend from the pythia expectation to no correlation 
between the decay electron and other hadrons. 



III. RESULTS & DISCUSSION 

Example jet functions, after efficiency corrections and 
combinatorial background subtraction, are shown in 
Fig. 2] for e,- nc , ebk g and enF triggers. Both the near- 
and away-side jet shapes are clearly present in the p+p 
data, however the statistical uncertainties in the Au+Au 
data are much larger. The boxes show the systematic un- 
certainties from all sources except for the overall normal- 
ization uncertainty of 7.9% in p+p and 9.4% in Au+Au. 
These jet functions, and others shown in Appendix IIVI 
are integrated and fit to extract the yields and widths 
that follow. Examples of the integrated near- and away- 
side conditional yields in p+p collisions are shown in 
Fig. [5] for the ej„ c , ebkg and enr triggers as a function 
of the hadron px ■ 



A. p+p Collisions 

The near-side (0< A0 < 1.25 rad) conditional yields of 
hadrons associated with heavy-flavor electrons are shown 
in Fig. [5] for the four electron px selections used in this 
analysis. While, in general, parton fragmentation favors 
the production of hadrons carrying a small fraction of 
the parent quark momentum, the heavy meson, D or 
B, resulting from a heavy quark fragmentation typically 
carries a large fraction of the heavy quark momentum 
(z = Ph ° dr ° n is peaked at «0.60 for charm and «0.85 
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FIG. 4: (color online) ei nc — h, e^kg — h and euF — h (solid circles) for p+p (top panel) and Au+Au (bottom panel) collisions 
for 2.0< pt,b <3.0 GeV/c and 1.5< pr,h <2.0 GeV/c. The overall normalization uncertain of 7.9% in p+p and 9.4% in Au+Au 
is not shown. 
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FIG. 5: (color online) Near (top) and away (bottom) side conditional yields in p+p collisions as a function of hadron pt- ei„ c 
triggers are shown as triangles, ebkg triggers are shown as open circles and eHP triggers are shown as solid circles. The boxes 
on the eHF — h points are the systematic uncertainties except for the overall normalization uncertainty of 7.9%, which is not 
shown. 
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for bottom) [45H49j . For comparison we also show the 
conditional yields from correlations between two charged 
hadrons from Ref. [14]. The spectra of near-side associ- 
ated hadrons is harder for the enF triggers than for the 
hadron trigg crs in all overlapping PT,trig 

selections. The 

near-side correlations are expected to be dominated by 
hadrons that are also from the decay of the heavy me- 
son. The large mass of the heavy meson translates to 
a wider expected near-side correlation when the hadron 
and the electron are both from the heavy meson decay. 
Figure [7] shows the Gaussian widths of the near-side con- 
ditional yields as a function of the associated hadron px- 
Also shown for comparison are the near-side widths for 
hadron-hadron correlations jl4| . which primarily come 
from light parton jets. The widths of the enr — h cor- 
relations are slightly wider for 2.0< enF <3.0 GeV/c, 
consistent with the near side being dominated by decay 
induced correlations. For higher px electrons the statisti- 
cal uncertainties become too large to make a quantitative 
statement. Results from powheg [28[, a next-to-leading- 
order Monte Carlo calculation, with charm and bottom 
contributions set by FONLL calculations [3l| are shown. 
These simulations are consistent with the data except for 
the lowest electron and hadron momenta. 

The away-side (summed over 2.51< A<j> < 7r rad) con- 
ditional yields are shown in Fig. [8j The yields on the 
away side are dominated by the fragmentation and de- 
cay of particles in the opposing jet. As discussed above 
the opposing dijet does not have to contain a balancing 
heavy-flavor quark; here the yields are a mix of heavy and 
light parton jets. At a given pt trigger bin the heavy- 
flavor electron triggered away-side spectrum is harder 
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FIG. 6: (color online) Near-side eHF — h conditional yields in 
p+p collisions (solid points) as a function of the associated 
hadron pr- For comparison hadron-hadron conditional yields 
in p+p collisions from Ref. [lij] are also shown (the A<p range 
for the hadron-hadron yields is A<f> < 7r/3). The associated 
hadron px spectra are harder for eHF — h than hadron-hadron 
conditional yields at the same pr,tri g range (the highest pr 
trigger selection is for 4.0-4.5 GeV/c for the eHF triggers and 
4.0-5.0 GeV/c for the hadron triggers). The overall normal- 
ization uncertainty of 7.9% is not shown. 



than the light hadron triggered one. In order to quantify 
the slope differences between light hadron triggers and 
heavy-flavor electron triggers we plot the inverse slope 
from an exponential fit of the away-side spectra in Fig. [5] 
However, the electron only carries a fraction of the 
heavy meson px- pythia was used to estimate the parent 
meson average px for both charm and bottom mesons and 
the results are shown in Table [TTTT When comparing the 
inverse slopes at similar meson px, as opposed to similar 
trigger particle px, the difference in the inverse slopes 
between the two trigger types largely disappears. 



B. Au+Au Collisions 

The conditional yields for the near-side (0< Acj) <1.25 
rad) for Au+Au collisions are compared to p+p collisions 
in Fig. rTU] for Au+Au collisions with 0-60% centrality. 
In addition, to provide a more direct comparison of the 
Au+Au and p+p conditional yields we construct the ratio 
of the conditional yield in Au+Au to p+p: 
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FIG. 7: (color online) Gaussian widths of near-side condi- 
tional yields as a function of the associated hadron pr for 
four pT, e selections. Solid circles show results from eHF — h 
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tions. Crosses are from powheg [2|| with charm and bottom 
combined according to FONLL calculations [3l| . 
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TABLE III: Mean transverse momentum of the parent D and B mesons contributing to the heavy-flavor electron pr bins used 
here. They are combined according to the fraction of heavy-flavor electrons from b quarks, ^ e _ > ^^_ >e ^ according to the FONLL 
calculations [3l[ (as shown in Ref. 29]) to determine the mean heavy meson transverse momentum. 
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FIG. 8: (color online) Away-side eHP — h (2.51< < ir) 
conditional yields in p+p collisions (solid points) as a func- 
tion of the associated hadron pr- For comparison hadron- 
hadron conditional yields in p+p collisions from Ref. [141 ] are 
also shown (the A(f> range for the hadron-hadron yields is 
2.51A<^) < 7r). The associated hadron pr spectra are harder 
for eHF — h than hadron-hadron conditional yields at the 
same pT,tri g range (the highest pr trigger selection is for 4.0- 
4.5 GeV/c for the eHF triggers and 4.0-5.0 GeV/c for the 
hadron triggers). The overall normalization uncertainty of 
7.9% is not shown. 



shown in Fig. [TT] In the absence of any nuclear ef- 
fects /aa will be unity. The near-side Jaa is consistent 
with one (x 2 /do f '=12.3/7, statistical uncertainties only). 
Naively, this might be expected since the near-side corre- 
lations in both p+p and Au+Au collisions are expected 
to largely be from the heavy meson decay. Since the 
decay length is long compared to the size and lifetime 
of the matter produced in Au+Au collisions the subse- 
quent decay of the heavy meson should be unmodified 
by the matter. However, it is possible that the charm 
and bottom contributions are altered from p+p collisions 
due to medium effects (such as different energy loss for 
charm and bottom quarks). Additionally, the measured 
hadrons are not solely from D and B decay, but also 
from the fragmentation of the heavy quarks and, possi- 
bly, from interactions between the heavy quark and the 
matter. Rather than attempt to disentangle these con- 
tributions (which would be highly model dependent) we 
leave it to theoretical models to reproduce the /aa with 
the combined hadron sources. 



FIG. 9: (color online) Inverse slope of the 6hf triggered away- 
side conditional hadron pr distributions shown in Fig. [8] as a 
function of pr of the trigger particle. For comparison fits to 
the hadron-hadron data from Ref. [1411 are also shown. 



We compare the enF — h Iaa values to those from 
hadron-hadron collisions at approximately the same me- 
son pt (see Table IIII[) . In hadron-hadron correlations 
the observed Jaa has a strong dependence on the pr 
of the trigger hadron [l4[. For 2< pr <3 GeV/c elec- 
trons the closest hadron px selection from Ref. [14[ was 
4< px <5 GeV/c and for 3< pr <4 GeV/c electrons it 
was 5< pr <10 GeV/c (see Table Hn). We observe the 
near-side 7aa for heavy-flavor electron triggers to be con- 
sistent with those from the comparison hadron triggered 
results, though the present uncertainties are too large 
to be sensitive to the excess seen in the hadron-hadron 
correlations. 

In order to be sensitive to possible modifications of 
the away-side jet shape, we measure the away-side yields 
in two A<fi ranges as shown in Fig. [12] for enF triggers 
with 2.0< p T ,e <3.0 GeV/c and 3.0< p T , e <4.0 GeV/c. 
The wide away-side range, 1.25< A(f> < tt rad is sensitive 
to the entire modified away-side shape and the smaller 
away-side range, 2.51< Aip < it rad, is sensitive to only 
the p+p like part of the away-side correlations. The ratio 
of conditional yields in Au+Au to p+p for both away- 
side A(j> ranges is shown in Fig. [T2] Iaa is largest at 
low hadron pt and decreases with increasing hadron pt- 
For comparison, the hadron-hadron Iaa from Ref. [14| is 
shown for the most closely matched meson pt selections 
(see Table IIIip . The heavy- flavor electron triggered Iaa 
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FIG. 10: (color online) Near-side (0< A(f> <1.25 rad) in- 
tegrated yield for Au+Au (solid circles) and p+p collisions 
(open circles) for 2.0< pt,e <3.0 GeV/c (top panel) and 
3.0< pr,e <4.0 GeV/c (bottom panel) as a function of the 
associated hadron pr- The overall normalization uncertainty 
of 9.4% in Au+Au and 7.9% in p+p is not shown. Points are 
slightly shifted horizontally for clarity. 



FIG. 11: (color online) Near-side (0< A(f> <1.25 rad) Iaa for 
2.0< p T ,e <3.0 GeV/c (top panel) and 3.0< p T , e <4.0 GeV/c 
(bottom panel) as a function of the associated hadron pr for 
6hf (solid points) and hadron (open points) triggers (from 
Ref. [1J). The gray band around unity shows the overall 
normalization uncertainty (12.4%), which moves all points 
together. Points are slightly shifted horizontally for clarity. 



is consistent with the hadron-hadron Iaa when compared 
at similar meson pt selections. 

This similarity could be due to the expected large frac- 
tion of gluons in the away-side distributions. If the away- 
side parton path lengths through the matter are similar 
between heavy-flavor electron and light hadron triggers 
then the corresponding Iaa values should also be simi- 
lar. However, some fraction of the away-side correlations 
should be due to correlated heavy-flavor quarks. This 
can be isolated in future measurements by triggering on 
back-to-back heavy-flavor electrons. We conclude that 
the present measurements are not sensitive to any differ- 
ences caused by back-to-back heavy-flavor pairs, either 
because the differences between away-side heavy-flavor 
electrons and light partons are small or because there 
are too few of them to significantly alter the Iaa values. 

Motivated by hadron-hadron correlations, we com- 
pared the away-side jet shape between p+p and Au+Au 
collisions. To quantify the shape differences we construct 
Rhs IM-i which is the yield/radian in the head region 



where the p+p jet is peaked (here 2.51 < A<j) < ir rad) 
divided by the yield/radian in the shoulder region where 
the enhancement in the Au+Au jet yield is observed in 
hadron-hadron correlations (1.25< A<fi <2.51 rad). The 
systematic uncertainties on the ratio are largely corre- 
lated between the head and shoulder region except for 
the uncertainty due to «2 in Au+Au collisions, which 
is anti-correlated because of the shape of the azimuthal 
modulation from V2 ■ In p+p collisions this ratio is large 
since the yield in the head region is much larger than the 
yield in the shoulder region. In hadron-hadron correla- 
tions for Au+Au collisions Rhs is observed to be smaller 
than in p+p collisions because of the increased yield in 
the shoulder region [14|. Fig. Q3] shows Rhs for enr — h 
correlations for 2.0< px.e <3.0 GeV/c as a function of the 
Pt of the associated hadron. Rhs is smaller for Au+Au 
collisions than for p+p collisions indicating that a similar 
away-side shape modification takes place for enF triggers 
as for hadron triggers (the head and shoulder A(f> regions 
are slightly different between this analysis and Ref. [l4| 
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FIG. 12: (color online) Away-side conditional yields for wide (left) and narrow (right) away-side A(j> integration ranges for 
Au+Au (solid points) and p+p (open points). Top panels show 2.0< pr,e <3.0 GeV/c and bottom panels shown 3.0< 
Pr,e <4.0 GeV/c. Upper limits are for 90% confidence levels. The overall normalization uncertainty of 9.4% in Au+Au and 
7.9% in p+p are not shown. Points are slightly shifted horizontally for clarity. 
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FIG. 13: (color online) 7aa as determined from the away-side yields in Fig. 1121 Two A<f> ranges are shown: 1.25< A(f> < n rad 
(left panels) and 2.51< A<f) < n rad (right panels). The gray band around unity shows the overall normalization uncertainty of 
12.4%, which moves all points together. For comparison hadron-hadron Taa values from Ref. [l4| are also shown for trigger pr 
selections where the parent heavy meson has similar pt to the trigger light hadron (see Table UTT)) . Points are slightly shifted 
horizontally for clarity. The solid horizontal line is at and the dashed horizontal line is at 1. 
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FIG. 14: (color online) Ratio of the yield in the head re- 
gion per radian to that in the shoulder region per radian for 
Au+Au (black) and p+p (red). 



preventing a direct comparison). No pt,h dependence of 
Rhs is observed, however the statistical uncertainties are 
quite large. 



IV. CONCLUSION & OUTLOOK 

Studies of the yields of electrons from the decay of 
heavy- flavor mesons in Au+Au collisions have challenged 
the picture of medium-induced radiative energy loss as 
the dominant mechanism by which high-py hadrons are 
suppressed. Correlations of hadrons from light quark and 
gluon jets have shown large modifications of the correla- 
tion patterns between p+p and Au+Au collisions. Study- 
ing the correlations of electrons from the decay of heavy 
mesons with other hadrons in the event provides more in- 
formation about how the charm and bottom quarks prop- 
agate through the matter and how the modified corre- 
lation structures observed in hadron-hadron correlations 
are produced. Thus, they are a crucial component of hard 
physics in relativistic heavy-ion collisions. The interpre- 
tation is complicated by the ambiguity in the away-side 
flavor and because the electron does not carry all of the 
parent meson's momentum. This makes understanding 
p+p collisions as a baseline very important. 

We have presented first measurements of the azimuthal 
correlations of electrons from heavy-flavor decay with 
hadrons in both p+p and Au+Au collisions. These mea- 
surements provide a first step in understanding correla- 
tions involving open heavy flavor in the hot matter pro- 
duced in heavy-ion collisions. The Gaussian widths of 
these correlations are consistent with expectations from 
simulations of the fragmentation of heavy quarks and 
the decay of heavy-flavor mesons, except for the low- 
est pt electrons and hadrons where some differences are 
observed. In p+p collisions the spectra of associated 
hadrons on both the near and away side are harder than 
in hadron-hadron correlations measured in the same trig- 



ger pt range. However, the level of away-side suppression 
at large px is consistent between electron and hadron 
triggers when the trigger charged hadron and the par- 
ent heavy meson are at approximately the same px ■ The 
ratio of yields in the head region to those in the shoul- 
der region decreases from p+p to Au+Au collisions in a 
manner qualitatively consistent with hadron-hadron col- 
lisions 14]. Further measurements sensitive to the par- 
tonic content of the away-side jets (heavy quarks or light 
quarks and gluons) are necessary to determine if this is 
due primarily to cases where the away-side parton is a 
light quark or gluon or if the suppression of away-side 
heavy-jet fragmentation is similar to those of light par- 
tons. 

Near future measurements of heavy- flavor triggered az- 
imuthal correlations hold particular promise. Data taken 
in 2010 has improved statistics and the HBD was success- 
fully operated allowing the rejection of some of the Dalitz 
and conversion electron background. Additionally, d+Au 
data taken by PHENIX in 2008 will help constrain any 
cold nuclear matter effects. Such effects are expected to 
be small at midrapidity, but are not well constrained by 
existing data. In future data taking, the silicon vertex 
detector will be installed, which will enable the separa- 
tion of electrons from D and B decay and increase accep- 
tance for measuring charged hadrons. Application of the 
techniques developed here on data taken with these up- 
grades in place will allow for more detailed heavy-flavor 
correlation measurements. 
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APPENDIX: Jet Functions 



Figures [15] and [16] show comparisons of euF — h jet 
functions for Au+Au and p+p collisions for the indicated 
electron triggers and hadron-px bins. Figures fTTH^IH show 
the ejjF — h jet functions for p+p collisions only for the 
indicated electron triggers and hadron-py bins. 
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FIG. 15: (color online) 6hf — h jet functions for Au+Au (solid blue circles) and p+p collisions for 2.0-3.0 GeV/c electron 
triggers and the hadron-pr bins indicated. 
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FIG. 16: (color online) 6hf — h jet functions for Au+Au (solid blue circles) and p+p collisions for 3.0-4.0 GeV/c electron 
triggers and the hadron-pr bins indicated. 




FIG. 17: (color online) chf — h jet functions in p+p collisions for 1.5-2.0 GeV/c electron triggers and the hadron-pr bins 
indicated. 
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FIG. 18: (color online) eHF — h jet functions in p+p collisions for 2.0-3.0 GeV/c electron triggers and the hadron-pT bins 
indicated. 
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FIG. 19: (color online) eHF — h jet functions in p+p collisions for 3.0-4.0 GeV/c electron triggers and the hadron-pr bins 
indicated. 
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FIG. 20: (color online) 6hf — h jet functions in p+p collisions for 4.0-4.5 GeV/c electron triggers and the hadron-pr bins 
indicated. 
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